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Abstract 

In this paper, we discuss an interaction between complex geome- 
try and integrable systems. Section 1 reviews the classical results on 
integrable systems. New examples of integrable systems, which have 
been discovered, are based on the Lax representation of the equations 
of motion. These systems can be realized as straight line motions on 
a Jacobi variety of a so-called spectral curve. In section 2, we study 
a Lie algebra theoretical method leading to integrable systems and we 
apply the method to several problems. In section 3, we discuss the 
concept of the algebraic complete integrability (a.c.i.) of hamiltonian 
systems. Algebraic integrability means that the system is completely 
integrable in the sens of the phase space being folited by tori, which in 
addition are real parts of a complex algebraic tori (abelian varieties). 
The method is devoted to illustrate how to decide about the a.c.i. of 
hamiltonian systems and is applied to some examples. Finally, in sec- 
tion 4 we study an a.c.i. in the generalized sense which appears as 
covering of a.c.i. system. The manifold invariant by the complex flow 
is covering of abelian variety. 

Mathematics Subject Classification (2000). 37.135, 70H06, 14H40, 14H70, 
14M10. 

1 Integrable systems 

Let M be an even-dimensional differentiable manifold. A symplectic struc- 
ture (or symplectic form) on M is a closed non-degenerate differential 2-form 
uj denned everywhere on M. The non-degeneracy condition means that 

Vx e M, V£ + 0, 3n : uj (£, n) + 0, (£, rj G T X M) . 

The pair (M, uj) is called a symplectic manifold. 

Example 1.1 The cotangent bundle T*M possesses in a natural way a sym- 
plectic structure. In a local coordinate (x\, ■ ■ ■ , x n , yi, . . . , y n ) , 2n = dim M, 
the form uj is given by uj = ^fc=i ^ x k A dyj,. 
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Example 1.2 Another important class of symplectic manifolds consists of 
the coadjoints orbits O C Q* , where Q is the algebra of a Lie group Q and 
G^ = {Ad* /J, : g G G} is the orbit of /j G G* under the coadjoint representa- 
tion. 

Theorem 1 a) Let I : T*M — ► T X M, u\ i — ► f, 6e a map defined by 
t^liv) = ^ (r},{;) , \/r] £ r x M. T/ien / is an isomorphism generated by the 
symplectic form lo. 

b) The symplectic form lo induces a hamiltonian vector field IdH : M — ► 
T X M, x i — ► IdH (x) , where H : M — > R, is a differ entiable function 
(called hamiltonian). In others words, the differential system defined by 

x(t) = X H (x (t)) = IdH (x) , 

is a hamiltonian vector field associated to the function H. The matrix that is 
associated to an hamiltonian system determine a symplectic structure. 

Proof, a) Denote by I~ l the map I~ l : T X M — ► T*M, £ i — ► I -1 (£) = lo^, 
with I" 1 (£) (77) = w| (77) = lo (77, £) , V77 G T X M. The fact that the form w is 
bilinear implies that 

+ = w fa, £1 + 6), 

= w (77, £1) +w (77,62) , 

= J -1 (£1) fa) + J -1 (6) fa), Vt/GT.M. 

Now, since dim T^M = dim T*M, to show that I" 1 is bijective, it suffices to 
show that is injective. The form lo is non-degenerate, it follows that 

Kerr 1 = {£ G T X M : lo (77, = 0, V77 G T X M} = {0} . 

Hence I -1 is an isomorphism and consequently I is also an isomorphism (the 
inverse of an isomorphism is an isomorphism) . 

b) Let (xi, . . . ,x m ) be a local coordinate system on M, (m = dimM). We 
have 

k=l K k=l K 

where I (dx^) = £ k G T X M is defined such that :V?7 G T X M, rj^ = dxj, (77) = 
lo (77, , (k-th component of^). Define (771, ... , rj m ) and . . . , to be 
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respectively the components of r\ and then 



(771, ...,ry m ) J -1 



V P 



where J 1 is the matrix defined by J 1 = [u ( 



. Since this 



matrix is invertiblqj, we can search £ fc such that : 

/ \ 



r 





1 ~k-th place 




V 

The matrix J -1 is invertible, which implies 



J 



Ve; 





f 





V°7 



from which £ fc =(k-th column of J), i.e., = Jjfc, 1 < i < m, and con- 
sequently £ fc = Y2iL 1 Jik-£^- It is easily verified that the matrix J is skew- 



indeed, it suffices to show that the matrix J 1 has maximal rank. Suppose this 
were not possible, i.e., we assume that ranfc(J _1 ) / m. Hence YHiLi fa§~; el - ) = 0' 
VI < j < m, with ai not all null and lo ^X^i ai a§~> a§~) — 0, VI < j < m. In fact, since 

uj is non-degenerate, we have YliLi a ' e§~ ~ 0- Now • • • , g§^) 1S a basis of T X M, 

then ai = 0, Vi, contradiction. 
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symmetricl. From (1) we deduce that 



h dx *h dx > ~h \h dx « dx *' 

k=l i=l i=l \k=l / 

Writing x(t) = YZi ffe « is seen that x i (*) = ELi ^g, 1 < * < 
j < m, which can be written in more compact form x(t) = J (x) this 
is the hamiltonian vector field associated to the function H. This concludes 
the proof of the theorem. 

We define a Poisson bracket (or Poisson structure) on the space C°° as 

{, } : C°° (M) x C°° (M) — ► C°° (M) , (F, G) i — ► {F, G} , 

where {F, G} = (i u F (X G ) = X G F (u) = uj(X g ,X f ). This bracket is 
skew-symmetric {F, G} = — {G, F} , obeys the Leibniz rule {FG, H} = 
F {G, H} + G {F, H} , and satisfies the Jacobi identity 

{{H, F},G} + {{F, G},H} + {{G, H},F} = 0. 

When this Poisson structure is non-degenerate, we obtain the symplectic 
structure discussed above. 



Consider now M = 1" x 1" and let p £ M. By Darboux's theorem [3], 
there exists a local coordinate system(xi, . . . ,x n ,y±, . . . ,y n ) in a neighbour- 
hood of p such that 



^ydxidyi dyidxij 



Then X H = ( ™£ - %* °) , and X H F = {H, F} , \/F G C» (M) . 



=1 \dxi dyi dm dxi 

A nonconstant function F is called an integral (first integral or constant of 
motion) of Xp, if XhF = 0. In particular, H is integral. Two functions 
F and G are said to be in involution or to commute, if {F, G} = 0. The 
hamiltonian systems form a Lie algebra. 

We now give the following definition of the Poisson bracket : 

dF dG 



\dx' dx I %J dxi dxj ' 

l >3 



2 Indeed, since w is symmetric i.e., u> ^gjr, gf-J = — ^ (af - ! gf - ) > it follows that J 1 
is skew-symmetric. Then, I = J.J -1 = (J _1 ) T .J T = —J -1 . J, and consequently J T = J. 
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After some algebraic manipulation, we deduce that If 

£ + + Jm^-) = 0, VI < < 2n, 

V dx k dx k dx k ) 

then J satisfies the Jacobi identity. 

Consequently, we have a complete characterization of hamiltonian vector 
field 

f) J-f 

x(t) = X H (x(t)) = J—,xeM, (2) 

where H : M — > R, is a differentiable function (the hamiltonian) and 
J = J (x) is a skew-symmetric matrix, possibly depending on x £ M, for 
which the corresponding Poisson bracket satisfies the Jacobi identity : 

{{H, F},G} + {{F, G},H} + {{G, H},F} = 0, 

with {H, F} = <§f , Jf> = £. . j y g |g, the Poisson bracket. 

Example 1.3 /In important special case is when J = ( ® I. J , where I 



I O 

is the nx n identity matrix. The condition on J is trivially satisfied. Indeed, 
here the matrix J do not depend on the variable x and we have 



1 j—1 i—i V ^^n^i dxi dxi dx n -\-i 

Moreover, equations (10) are transformed into 

. _ dH . _ dH . dH . dH 

dpi dp n dqi dq n 

oil qi = xi, . . . , q n = x n ,p\ = x n+ i, . . . ,p n = xm- These are exactly the well 
known differential equations of classical mechanics in canonical form. 

It is a fundamental and important problem to investigate the integrability 
of hamiltonian systems. Recently there has been much effort given for finding 
integrable hamiltonian systems, not only because they have been on the 
subject of powerful and beautiful theories of mathematics, but also because 
the concepts of integrability have been applied to an increasing number of 
applied sciences. The so-called Arnold-Liouville theorem play a crucial role 
in the study of such systems; the regular compact level manifolds defined by 
the intersection of the constants of motion are diffeomorphic to a real torus 
on which the motion is quasi-periodic as a consequence of the following purely 
differential geometric fact : a compact and connected n-dimensional manifold 
on which there exist n vector fields which commute and are independent at 
every point is diffeomorphic to an n-dimensional real torus and each vector 
field will define a linear flow there. 
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Theorem 2 (Arnold-Liouville theorem)[3, 14] : Let Hi = H, H2, H n , be 
n first integrals on a 2n- dimensional symplectic manifold that are functionally 
independent (i.e., dH\ A ... A dH n 7^ 0), and pairwise in involution. For 
generic c = (c±, c n ) the level set 

n 

M c = p| {x e M : Hi (x) = a, a € M} , 

i=l 

will be an n-manifold. If M c is compact and connected, it is diffeomorphic to 
an n-dimensional torus T n = R n /Z n and the solutions of the system (2) are 
then straight-line motions on T n . If M c is not compact but the flow of each 
of the vector fields Xn k is complete on M Cl then M c is diffeomorphic to a 
cylinder M fc x T n ~ k under which the vector fields Xjj k are mapped to linear 
vector fields. 

As a consequence, we obtain the concept of complete integrability of a 
hamiltonian system. For the sake of clarity, we shall distinguish two cases : 

a) Case 1 : det J / 0. The rank of the matrix J is even, m = 2n. A hamilto- 
nian system (2) is completely integrable or Liouville-integrable if there exist 
n firsts integrals H\ = H,H2, ■ ■ ■ ,H n in involution, i.e., {Hk,Hi} = 0, 1 < 
k,l < n, with linearly independent gradients, i.e., dH\ A ... A dH n 7^ 0. For 
generic c = (c\, c n ) the level set 

n 

M c =f]{xeM:Hi{x)= c u Ci G R} , 
i=i 

will be an n-manifold. By the Arnold-Liouville theorem, if M c is compact 
and connected, it is diffeomorphic to an n-dimensional torus T n = M n /Z n 
and each vector field will define a linear flow there. In some open neigh- 
bourhood of the torus there are coordinates s±, . . . , s n , ip\, . . . , (p n in which 
u takes the form u> = J2k=i d s k A difk- Here the functions s^ (called action- 
variables) give coordinates in the direction transverse to the torus and can 
be expressed functionally in terms of the firsts integrals H^. The functions 
(fk (called angle-variables) give standard angular coordinates on the torus, 
and every vector field Xn k can be written in the form tpk = (si, . . . , s n ) , 
that is, its integral trajectories define a conditionally-periodic motion on the 
torus. In a neighbourhood of the torus the hamiltonian vector field Xn k 
take the following form = 0, cp^ = (si, . . . , s n ) , and can be solved by 
quadratures. 

b) Case 2 : det J = 0. We reduce the problem to m = 2n + k and we look for 

k Casimir functions (or trivial invariants) H n+ \, H n+ k, leading to identi- 

OH 

cally zero hamiltonian vector fields J =0, 1 < i < k. In other words, 
the system is hamiltonian on a generic symplectic manifold 

n+k 

f] {x G M. m : Hi (x) = C{} , 
i=n+l 
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of dimension m — k = 2n. If for most values of Cj 6 R, the invariant manifolds 

n+fc 

P| g R m : tf, (x) = Ci} , 
i=i 

are compact and connected, then they are n-dimensional tori T n = R n /Z n 
by the Arnold-Liouville theorem and the hamiltonian flow is linear in angular 
coordinates of the torus. 

2 Isospectral deformation method 

A Lax equation is given by a differential equation of the form 

A(t) = [A(t),B(t)} or [B(t),A(t)}, (3) 

where 

N N 

A(t) = J2 A k (*) h k , B(t) = J2 B k (t) h\ 

k=l k=l 

are functions depending on a parameter h (spectral parameter) whose coeffi- 
cients Ak and are matrices in Lie algebras. The pair (A, B) is called Lax 
pair. This equation established a link between the Lie group theoretical and 
the algebraic geometric approaches to complete integr ability. The solution 
to (3) has the form A(t) = g(t)A(0)^(f) _1 , where g(t) is a matrix defined as 
g (t) = —A(t)g(t). We form the polynomial P (h, z) = det (A — zl) , where 
z is another variable and I the n x n identity matrix. We define the curve 
(spectral curve) C, to be the normalization of the complete algebraic curve 
whose affine equation is P (h, z) = 0. 

Theorem 3 The polynomial P (h, z) is independent of t. Moreover, the func- 
tions tr (A n ) are first integrals for (3). 

Proof. Let us call L = A — zl. Observe that 

P = det L.tr (L^L) = det L.tr {L~ X BL - B) = 0, 
since trL~ l BL = trB. On the other hand 

An = AA 11 - 1 + AAA n ' 2 + --- + A n ~ 1 A, 

= [A, B] A n - 1 + A [A, B] A n ~ 2 + ■■■ + A 11 " 1 [A, B] , 

= (AB - BA) A" 1 ' 1 + ■■■ + A n - 1 (AB - BA) , 

= ABA 11 ' 1 - BA n + ■■■ + A n B - A^BA, 

= A (BA" 1 - 1 ) - (BA" 1 - 1 ) A + ■ ■ ■ + A (A n_1 B) - (A^B) A. 
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Since tr (X + Y) = trX + trY, trXY = trYX, X,Y G M n (C) , we obtain 



and consequently tr (A n ) are first integrals of motion. This ends the proof 
of the theorem. 



We have shown that a hamiltonian flow of the type (3) preserves the 
spectrum of A and therefore its characteristic polynomial. The curve C : 
P(z,h) = det (A(h) — zT) = 0, is time independent, i.e., its coefficients 
tr (A n ) are integrals of the motion (equivalently, A(t) undergoes an isospec- 
tral deformation. Some hamiltonian flows on Kostant-Kirillov coadjoint or- 
bits in subalgebras of infinite dimensional Lie algebras (Kac-Moody Lie al- 
gebras) yield large classes of extended Lax pairs (3) . A general statement 
leading to such situations is given by the Adler-Kostant-Symes theorem. 

Theorem 4 Let £ be a Lie algebra paired with itself via a nondegenerate, 
ad-invariant bilinear form ( , ), C having a vector space decomposition C = 
K, + Af with K, and Af Lie subalgebras. Then, with respect to ( ,), we have 
the splitting C = C* = K, L + Af 1 - and Af* = K, 1 - paired with Af via an induced 
form ((,)) inherits the coadjoint symplectic structure of Kostant and Kirillov; 
its Poisson bracket between functions H\ and H2 on Af* reads 

{H 1 ,H 2 }(a) = ((a,[V^H 1 ,V^H 2 })} , a G M* . 

Let V C M* be an invariant manifold under the above co-adjoint action of 
Af on Af* and let AiV) be the algebra of functions defined on a neighborhood 
of V , invariant under the coadjoint action of C (which is distinct from the 
Af — Af* action). Then the functions H in AiV) lead to commuting Hamil- 
tonian vector fields of the Lax isospectral form 

a = [a,prtc(VH)] , pr^ projection onto K. 

This theorem produces hamiltonian systems having many commuting inte- 
grals ; some precise results are known for interesting classes of orbits in both 
the case of finite and infinite dimensional Lie algebras. Any finite dimensional 
Lie algebra C with bracket [, ] and killing form (, ) leads to an infinite dimen- 
sional formal Laurent series extension C = YZoo A i h% N £Z free, 
with bracket [Yl Aih 1 , ^ Bjh-*] = Ylij ^ i+J > and ad-invariant, sym- 
metric forms (^2 Aih 1 , Bjh-*) k = J2i+j=-k (Ai,Bj) , depending on k G Z. 
The forms (, } k are non degenerate if (, ) is so. Let (p < q) be the vector 
space of powers of h between p and q . A first interesting class of problems 
is obtained by taking C = Gl(n,M) and by putting the form on the 
Kac-Moody extension. Then we have the decomposition into Lie subalge- 
bras C = £ ,oo + £-00,-1 = K, + Af with K = K.- 1 , Af = Af 1 and K = Af*. 
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Consider the invariant manifold V m , m > 1 in /C = N* , denned as 

f m -1 1 

V m = < A = Aih 1 + ah m , a = diag(ai, • • • , a n ) fixed > , 
with diag (A m -i) = 0. 

Theorem 5 The manifold V m has a natural symplectic structure, the func- 
tions H = (f(Ah~i), on V m for good functions f lead to complete inte- 
grable commuting hamiltonian systems of the form 

m—l 

A = ]A,pr K {f'{Ah-i)h k ->^ , A = ^ A i h i + ah, 

i=0 

and their trajectories are straight line motions on the jacobian of the curve 
C of genus (n - 1) (nm - 2) /2 defined by P (z, h) = det (A - zl) = 0. The 
coefficients of this polynomial provide the orbit invariants of V m and an in- 
dependent set of integrals of the motion ( of particular interest are the flows 
where j = m, k = m + 1 which have the following form, 

A = [A , adp ad^Am-x + @h] , A = /' (ccj) , 

the flow depends on f through the relation Pi = f (a^) only). 

Another class is obtained by choosing any semi-simple Lie algebra L . Then 
the Kac-Moody extension C equipped with the form (, ) = (, ) has the nat- 
ural level decomposition C = Yliez [Li,Lj] C [Lo,A>] = 0> L* = 
L_j. Let B + = J2i>o^i an d B~ = Yltio L-i ■ Then the product Lie algebra 
C x C has the following bracket and pairing 

\h,h),(l'iX)] = (m'i],-M;>]), ((h,h),(U2)) = (h,t 1 )-(hX)- 

It admits the decomposition into K, + N with 

K = {(l,-l):l€C}, K L = {(/,/):/ G £} , 

JV={(Z_,Z+) :/_ e£- l + €B + ,jw (i-)=Pr , oG+)}, 

= {(/_,/+) :/_ GB-,Z + G5+,pr (Z + + /-)=0}, 

where pro denotes projection onto Lq. Then from the last theorem , the 
orbits in J\f*=K ± possesses a lot of commuting hamiltonian vector fields of 
Lax form: 

Theorem 6 The N-invariant manifolds V-j,k = Yl-j<i<k Q £ ~ /C^ - , 
has a natural symplectic structure and the functions H{l\,l2) = f{h) on 
V-j t k lead to commuting vector fields of the Lax form 



I 



I, {pr + - -pr )VH 



pr + projection ontoB' 



their trajectories are straight line motions on the Jacobian of a curve defined 
by the characteristic polynomial of elements in V-j t k. 



9 



Using the van Moerbeke-Mumford linearization method [21], Adler and 
van Moerbeke [1] showed that the linearized flow could be realized on the 
jacobian variety Jac(C) (or some sub-abelian variety of it) of the algebraic 
curve (spectral curve) C associated to (3). We then construct an algebraic 
map from the complex invariant manifolds of these hamiltonian systems to 
the jacobian variety Jac{C) of the curve C. Therefore all the complex flows 
generated by the constants of the motion are straight line motions on these 
jacobian varieties i.e. the linearizing equations are given by 



where . . ,co g span the g- dimensional space of holomorphic differentials 
on the curve C of genus g. In an unifying approach, Griffiths [8] has found 
necessary and sufficient conditions on B for the Lax flow (3) to be linearizable 
on the jacobi variety of its spectral curve, without reference to Kac-Moody 
Lie algebras. 

Next I schall discuss a number of integrable hamiltonian systems. 
2.1 The Euler rigid body motion. 

It express the free motion of a rigid body around a fixed point. Let M = 
(mi, m2, m 3 ) be the angular momentum, Q = (mi /Ii,m,2/l2, mz/Iz) the 
angular velocity and ii, 1% et ^3, the principal moments of inertia about the 
principal axes of inertia. Then the motion of the body is governed by 



If one identifies vectors in M 3 with skew-symmetric matrices by the rule 




M=MAfi. 



(4) 




-a 3 




-ai 





then a A b 1 — ► [A, B] = AB — BA. Using this isomorph 
and (so(3), [, ]), we write (4) as M = [M,Q] , where 



.ism between (M 3 , A) 




Now m = m 



this implies that 



M = [M, AM] , 



(5) 
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where 

AM = I A 3 m 3 -Aimi ) Gso(3), 
\ -A 2 m 2 Aimi 

with Aj = If 1 . Equation (5) is explicitly given by 

rn\ = (A 3 - A 2 ) m 2 m, 3 , 

rh 2 = (Ai-A 3 )mim 3 , (6) 
= (A 2 -Ai)mim 2 , 

and can be written as a hamiltonian vector field 

x = J—, x = (mi,m 2) m 3 ) 1 , 
with the hamiltonian H = ^ (Aimf + A 2 m,2 + A3m|) , and 

— 771-3 m 2 

J = | m 3 —mi ) 6 so (3). 
-m 2 mi 

We have det J = 0, so m = 2n + and m — k = rk J. Here m = 3 and rfe 
J = 2, then n = fe = 1. The system (6) has beside the energy Hi = H, a 
trivial invariant H2, i.e., such that: «/^jjr = 0, or 

/ Mi 

dm 1 
<9H 2 
9m2 
9g 2 
\ <9m 3 

implying = mi, = m 2 , = m 3 , and consequently 

The system evolves on the intersection of the sphere Hi = ci and the el- 
lipsoid H2 = c 2 . In M 3 , this intersection will be isomorphic to two circles 
^with || < ci < . We shall show that the problem can be integrated in 
terms of elliptic functions, as Euler discovered using his then newly invented 
theory of elliptic integrals. Observe that the first equation of (6) reads 

dm-\ , . , 

(A 3 - A 2 ) dt, (7) 





m 2 m 3 

where mi,m 2 and WI3 are related by 



Aim-i + A 2 m| + AsmJ = ci, mf + m| + m 3 = c 2 . 
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Therefore, if A2 7^ A3, we have 



_ ± / c 2 A 3 - ci + (Ai - A 3 )m 2 m - ± Cl ~ C2 ^ 2 + ~ 
y A 3 — A2 y A3 — A2 

Substituting these expressions into (7), we find after integration that the 
system (6) amounts to an elliptic integral 

/•mi(t) dm 
I = = ct, 

imi(o) v (m 2 + a) (m 2 + 6) 

with respect to the elliptic curve 

C : w 2 = (z 2 + a) (z 2 + b) , (8) 

with a = c \* 3 _~xf , b = X 2 -t L 2 > c = - A 3 ) (A 2 - Ai). Then the func- 

tions rriiit) can be expressed in terms of theta-functions of t, according to 
the classical inversion of abelian integrals. 

We shall use the Lax representation of the equations of motion to show 
that the linearized Euler flow can be realized on an elliptic curve isomorphic 
to the original elliptic curve (8). The solution to (5) has the form 

M (t) = 0(t)M (t) M T (t) , 

where O (t) is one parameter sub-group of SO (3) . So the hamiltonian flow 
(5) preserves the spectrum of X and therefore its characteristic polynomial 
det (M — zl) = — z (z 2 + m\ + m 2 + m 3 ) . Unfortunately, the spectrum of a 
3x3 skew- symmetric matrix provides only one piece of information; the 
conservation of energy does not appear as part of the spectral information. 
Therefore one is let to considering another formulation. The basic obser- 
vation, due to Manakov [20] , is that equation (5) is equivalent to the Lax 
equation 

A = [A,B], 

where A = M + ah, B = AM + @h, with a formal indeterminate h and 




a 



Ai = 

Ctl — CK3 «2 — Oil 

and all distinct. The characteristic polynomial of A is 

P (h, z) = det (A - zl) , 

= det (M + ah - zl) , 

3 / 3 

= H( aj h-z)+ k>. 
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The spectrum of the matrix A = M + ah as a function of h € C is time 
independent and is given by the zeroes of the polynomial P (h, z) , thus 
defining an algebraic curve (spectral curve). Letting w = h/z, we obtain the 
following elliptic curve 

3 

z 2 Yl (ctjw - 1) + 2H lW - 2H 2 = 0, 

3=1 

which is shown to be isomorphic to the original elliptic curve. Finally, we 
have the 

Theorem 7 The Euler rigid body motion is a completely integrable system 
and the linearized flow can be realized on an elliptic curve. 

2.2 The geodesic flow for a left invariant metric on SO (4) . 

Consider the group £0(4) and its Lie algebra so(4) paired with itself, via 
the customary inner product (X,Y) = — ^ tr (X.Y) , where 



X 



-X 3 

x 3 

-X 2 Xl 



X2 
-Xl 





-X 4 \ 



-X5 





G so(4). 



A left invariant metric on SO (4) is defined by a non-singular symmetric 
linear map A : so(4) — ► so(4), X i — > AX, and by the following inner 
product; given two vectors gX and gY in the tangent space SO (4) at the 
point g £ SO (4), {gX,gY) = (X, A^.Y) . Then the geodesic flow for this 
metric takes the following commutator form (Euler- Arnold equations) : 



X = [X, AX] , 



(9) 



A.X 



where 

-A3X3 

A3X3 

-X2X2 X1X1 

\ A4X4 A5X5 

In view of the isomorphism between 
system (9) as 



A2Z2 
-A1X1 


A 6 x 6 



-A4X4 \ 
-A5X5 
-A 6 x 6 




G so(4). 

, a) , and (so (4) , [, ]) we write the 



xi = (A 3 - A 2 ) x 2 x 3 + (A 6 - A 5 ) x 5 x 6 , 

*2 = (Ai - A 3 ) X1X3 + (A 4 - A 4 ) x 4 x 6 , 

*3 = (A2 - Ai) xix 2 + (A 5 - A 4 ) X4X5, 

±4 = (A 3 - A 5 ) X3X5 + (A 6 - A 2 ) x 2 x 6 , 

x 5 = (A 4 - A 3 ) X3X4 + (Ai - A 6 ) xix 6 , 

x$ = (A 2 - A 4 ) x 2 x 4 + (A 5 - Ai) X1X5. 
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These equations can be written as a hamiltonian vector field 



with 



H 



T 9H 



{X, AX) = - (Aire? + \ 2 x\ + ■■■ + \qx 2 6 ) , 



(10) 



the hamiltonian and 



J 



( 





~X3 


X2 





-Xq 


X5 


\ 




X3 





-Xl 


Xq 





— X4 






-X2 


Xl 





-x 5 


X4 












-Xq 


X5 





— X3 


X2 






Xq 





— X4 


^3 





-xi 




V 


-%5 


X4 





~X2 


Xl 





/ 



G so(6). 



We have det J = 0, so m = 2n + and m — k = rk J. Here to = 6 and rg 
J = 4, then n = k = 2. The system (10) has beside the energy Hi = H, two 
trivial constants of motion : 



H 2 
H 3 



1 



>? + xl + • ■ ■ + x|) 



X1X4 + x 2 x 5 + X 3 X 6 . 



Recall that #2 and H3 are called trivial invariants (or Casimir functions) 
because </^r = J^fjr = 0- In order that the hamiltonian system (10) be 
completely integrable, it is suffices to have one more integral, which we take 
of the form 

#4 = 2 + w x l ^ — i" ^ x D ■ 

The four invariants must be functionally independent and in involution, so 
in particular 



i.e., 



((A 3 - A 2 ) m + (Ai - A 3 ) H2 + (A 2 - Ai) fi 3 ) xix 2 x 3 
+ ((A6 - A 5 ) m + (Ai - A 6 ) fi 5 + (A 5 - Ai) hq) xix 5 xq 

A2) M4 + (A 2 - A 4 ) hq) X2X4XQ 

A5) + (A4 - A3) fl 5 ) X3X4X5 = 0. 



+ ((A 4 - A 6 )^ 2 + (A 6 
+ ((A 5 -A 4 )^3 + (A3 



Then 



(A 3 - A 2 ) m + (Ai 
(A6 - A 5 ) m + (Ai 
(A 4 - A 6 ) jji 2 + (A 6 



As)/^2 + (A 2 - Ai)/x 3 = 0, 
A6)^5 + (A5 - Ai)/x 6 = 0, 
A2) M4 + (A 2 - A 4 ) hq = 0, 
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(A 5 - A 4 ) ^3 + (A 3 - A 5 ) /x 4 + (A 4 - A 3 ) /u 5 = 0. 



Put 



A 



( A 3 - A 2 Ai - A 3 A 2 - Ai \ 

A 6 - A 5 Ai - A 6 A 5 - Ai 

A 4 -A 6 A 6 -A 2 A 2 -A 4 

V A 5 - A 4 A 3 - A 5 A 4 - A 3 / 

The number of solutions of this system is equal to the number of columns 
of the matrix A minus the rank of A. If rkA = 4, we have two solutions : 
Hi = 1 lead to the invariant H 2 and [ii = Aj lead to the invariant H%. This 
is unacceptable. If rkA = 3, each four-order minor of A is singular. Now 

/ A3 — A 2 Ai — A3 A 2 — Ai \ 



A 6 - A 2 

V A 5 - A 4 A 3 - A 5 / 

\ 

Ai-A 6 
A 6 - A 2 

A3 — A5 A 4 — A3 / 

\ 

Ai — A6 A5 — Ai 
A 2 - A 4 



A 6 - A 5 
A 4 - A 6 


/ Ai — A3 A 2 — Ai 

A 4 - A 6 

A 5 - A 4 

/ A 2 - Ai 

A 6 - A 2 



V 



\ A5 — A 4 A3 — A5 A 4 — A3 



— — (A6 - A5) C, 



(Ai - A 6 ) C, 



(A 2 - Ai) C, 







J 



where 
C = 



AiAgA 4 + A1A2A5 — AiA 2 A 4 + A3A6A5 — A3AgA 4 — A3A2A5 
+A 4 A 2 A5 + A 4 AiA3 — A 4 AiA5 + AeA 2 A3 — AeA 2 As — A1A6A3, 



and it follows that the condition for which these minors are zero is C = 0. 
Notice that this relation holds by cycling the indices : Ia4,2rv5,r\3a 
6. Under Manakov [20] conditions, 



Ai = 
A 4 = 



Pi 


-03 


A 2 = 


Pi 


-P 3 


A - & 
A3 — 


-Pi 


Oil 


- 03' 




ai 


- «3 ' 




— a 2 


Pi 


-Pa 


A 5 = 


Pi 


-Pi 


A - A 
A 6 — 


~P4 


Ct\ 


— a 4 ' 




a 2 


— a 4 ' 


03 


— a 4 



(11) 



where ai,Pi £ C, Y\ i<C j (a« — Pj) / 0, equations (10) admits a Lax equation 



with an indeterminate h : 



a 



(X + ah) = [X + ah, AX + j3h] , 
( an ^ ( p! \ 

p 2 



a 2 
a 3 

\ Q 4 / 



P 



p 3 
\ p A J 



(12) 
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t 

X = [X, A.X] & (9) , 
[X,P] + [a, A.X] =0^ (11), 
[a, j3\ = trivially satisfied for diagonal matrices. 

The parameters //i, . . . , can be parameterized (like Ai, . . . , Xq) by : 

72 - 73 71 - 73 7i - 72 

Ml = , M2 = , M3 



OL2 — «3 ai — CK3 «i — Oi2 



7i - 74 72 - 74 73 - 74 

M4 = , = , M6 



Oi\ — 0!4 Ot2 — CK4 Us — «4 

To use the method of isospectral deformations, consider the Kac-Moody 
extension (n = 4): C = j^-oo^^ : A arbitrary G Z, G gZ(n,M)j , of 

^(n,R) with the bracket: [E^,E^] = Zk (E i+j =k [A, B;]) ^, 
and the ad-invariant form: Aih 1 ,^2 Bjh^) = E i+ , = _ 1 (ij,-Bj}, where 
(,) is the usual form defined on gZ(n,IR). Let /C and A/" be respectively the 
> and < powers of h in C, then £ = K. + M, for the pairing defined 
above JC = JC L , N = so that JC = AT*. The orbits described in 

this way come equipped with a symplectic structure with Poisson bracket 
{H 1 ,H 2 }(a) = {a,[V K *H l ,V K *H 2 }) , where a G JC* and V K *H G JC. Ac- 
cording to the Adler-Kostant-Symes theorem, the flow (12) is hamiltonian 
on an orbit through the point X + ah, X G so(4)) formed by the coadjoint 
action of the subgroup Gj\f C SL (n) of lower triangular matrices on the dual 
Kac-Moody algebra M* ~ /C 1 - = JC. As a consequence, the coefficients of 
z l h l appearing in curve : 

r : {(z,h) GC 2 : det (X + ah - zl) = 0} , (13) 

associated to the equation (12) , are invariant of the system in involution for 
the symplectic structure of this orbit. Notice that 

det (gXg -1 ) = det X = (xix^ + x 2 x^ + x^xq) 2 , 

tr {gXg- 1 ) 2 = tr (gX 2 g- 1 ) = tr (A 2 ) = -2 (x 2 + x\ + • • • + x 2 ) . 

Also the complex flows generated by these invariants can be realized as 
straight lines on the abelian variety defined by the periods of curve 17 Ex- 
plicitly, equation (13) looks as follows 

4 

T : Yl ( a i h - z ) + 2H ± h2 - 2H lZ h + 2H 2 z 2 + tff = 0, 



16 



where H^X) = ci, H 2 {X) = c 2 , H 3 (X) = 2H = c 3 , H 4 (X) = c 4 . with 
c i) c 2 5 C3,C4 generic constants. T is a curve of genus 3 and it has a natural 
involution a : F — > F , (z,h) t— > (—z,—h). Therefore the jacobian variety 
Jac(T) of T splits up into an even and old part : the even part is an el- 
liptic curve To = F/a and the odd part is a 2— dimensional abelian surface 
Prym{F /Fq) called the Prym variety : JaciF) = Fq + PrymiF /Fq). The 
van Moerbeke-Mumford linearization method provides then an algebraic map 
from the complex affine variety Hi=i {Hi{X) = q} C C 6 to the Jacobi vari- 
ety Jac(F). By the antisymmetry of F , this map sends this variety to the 
Prym variety Prym (F/Fq) : 

4 3 
p| {Hi(X) = a} -> Prym(F /F ) , p -> £ s k , 

i=l k=l 

and the complex flows generated by the constants of the motion are straight 
lines on Prym(F /Fq). Finally, we have the 

Theorem 8 The geodesic flow (9) is a hamiltonian system with 

h = h x = x - {\ x x\ + \ 2 x\ + ■■■ + Xexj) , 
the hamiltonian. It has two trivial invariants 

H 2 = -{a$ + a$ + --- + 4), 

H 3 = xix 4 + x 2 x 5 + x 3 x 6 . 

Moreover, if 

A1A6A4 + A1A2A5 — A1A2A4 + A3A6A5 — A3A6A4 — A3A2A5 
+A4A2A5 + A4A1A3 — A4A1A5 + A6A2A3 — A6A2A5 — A1A6A3 = 0, 

the system (9) has a fourth independent constant of the motion of the form 

H ± = \ + ^2^2 H h l^xf) . 

Then the system (9) is completely integrable and can be linearized on the 
Prym variety Prym(F /Fq). 

2.3 The Toda lattice. 

The Toda lattice equations (discretized version of the Korteweg-de Vries 
equatiorJl) motion of n particles with exponential restoring forces are gov- 
erned by the following hamiltonian 

N N 

1=1 1=1 

3 In short K-dV equation : |^ — 6u|^ + |pf = 0. This is an infinite-dimensional 
completely integrable system. 
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The hamiltonian equations can be written as follows 



Pi, Pi 



In term of the Flaschka's variables [6] : 
equations take the following form 



IpQi-Qi+l h. 



Oj (b i+ i - bi) , hi = 2(o, 



a l-i), 



\pi, Toda's 



(14) 



with b]\r+i = b\ and ao = ajy. To show that the system (14) is completely 
integrable, one should find N first integrals independent and in involution 
each other. Prom the second equation, we have 



i N N „ 

d v-^ s-^dbi 



i=i 



dt 



0, 



and we normalize the 6j's by requiring that = 0- This is a first 

integral for the system. We further define N x N matrices A and B with 



A 



( h 

cti 




\ a N 



ai 
b 2 







a N 







bw-i cln-i 
a-N-i bN ) 



,B = 



( o 

\ a N 



ai 




-a N \ 



—a-N-i 







Then (14) is equivalent to the Lax equation 

A = [B,A}. 

From theorem 3, we know that the quantities Ik = \trA k , 1 < k < N, are 
first integrals of motion : To be more precise 

I k = tr{A.A k ~ 1 ) = tr([B,A].A k - 1 ) = tr(BA k - ABA^ 1 ) = 0. 

Notice that I\ is the first integral already know. Since these N first integrals 
are shown to be independent and in involution each other, the system (14) 
is thus completely integrable. 



2.4 The Garnier potential. 

Consider the hamiltonian 



H = \ (*? + xl) - \ (X iy f + X 2 yl) + \(yl + ytf , (15) 



where Ai and A2 are constants. The corresponding system is given by 

m = xi, x x = (Ai -yl-yl)yi, (16) 
m = x 2 , x 2 = (A 2 



vl) yi, 
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Theorem 9 The system (16) has the additional first integral 

H 2 = j (( x i^2 - x 2 yi) 2 - (X 2 yf + Aiy|) - (Ai + A 2 ) y\yfj 
+\ (AiA 2 (yj + yl) - {X 2 x\ + Ai^)) . 

and is completely integrable. The flows generated by Hi = H(15) and H 2 
are straight line motions on the jacobian variety of a smooth genus two hy- 
perelliptic curve 71(17) associated to a Lax equation. 

Proof. We consider the Lax representation in the form A = [A, B] , with the 
following ansatz for the Lax operator 

where 

V - _( ft -A 1 )( ft -A 2 )(l + I(^ + ^)), 

U = 2 (/l " Al)(/l - A2)( /^ + /^ ) ' 

W = (/ l - Al )(/ l -A 2 )(l(^ x - + ^-)- h +\{yl + yl)), 

R = h-yl-yl. 

We form the curve in (z,h) space 

P{h, z) = det (A - zl) = 0, 

whose coefficients are functions of the phase space. Explicitly, this equation 
looks as follows 

n-.z 2 = p 5 (h), (17) 

= (h - Ai)(/i - X 2 )(h 3 - (Ai + \ 2 )h 2 + (AiA 2 - H-fih - H 2 ), 
with H\ (15) the hamiltonian and a second quartic integral H 2 of the form 

H 2 = -j{X 2 yf + Xiy 2 + (Xi + \ 2 )y\y\ - (xiy 2 - x 2 yi) 2 ) 
~(X 2 x 2 + Xix 2 2 - XiX 2 (yj + y 2 2 )). 

The functions Hi and H 2 commute : {Hi,H 2 } = and the system (16) is 
completely integrable. The curve 7i determined by the fifth-order equation 
(17) is smooth, hyperelliptic and its genus is 2. Obviously, 7i is invariant 
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under the hyperelliptic involution (h, z) rx (h, —z) . Using the van Moerbeke- 
Mumford linearization method, we show that the linearized flow could be 
realized on the jacobian variety Jac (H) of the genus 2 curve TL. For generic 
c = (01,02) £ C 2 the afhne variety defined by 

2 

M c = P| {x g C 4 : H t (x) = a} , 
i=i 

is a smooth affine surface. According to the schema of [5] , we introduce 
coordinates s\ and s 2 on the surface M c , such that M c (sj) = 0, Ai / A2, i.e., 

si + s 2 = - (yj + yl) + Ai + A 2 , sis 2 = - (\ 2 yj + Xiyf) + AiA 2 . 

After some algebraic manipulations, we obtain the following equations for s\ 
and s 2 : 

Si — 2 , S 2 — ^ , 

si - s 2 s 2 - Si 

where P5 (s) is defined by (17) . These equations can be integrated by the 
abelian mapping 

n — ► jac(w) =c 2 /l , P ^ ( r ui , r^ 2 ), 

where the hyperelliptic curve TC of genus two is given by the equation (17), 
L is the lattice generated by the vectors n\ + J7n 2 ,(ni,ra 2 ) G Z 2 ,S1 is the 
matrix of period of the curve TL, (wi,u; 2 ) is a canonical basis of holomorphic 
differentials on H, i.e., 

ds sds 



and po is a fixed point. This concludes the proof of the theorem. 

2.5 The coupled nonlinear Schrodinger equations. 

The system of two coupled nonlinear Schrodinger equations is given by 

i ^ z +^ + n a+ 2 -(\a\ 2 + \b\ 2 )a+ 1 -(a* + bi)a = 0, (18) 

^ + S-^ + ^( H2 + |6|2 ) 6+ ^ a2 + 62 )' = °' 

where a(z,t) and b(z,t) are functions of z and t, the bar "— " denotes the 
complex conjugation, "||" denotes the modulus and is a constant. These 
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equations play a significant role in mathematics, with a important number 
of physical applications. We seek solutions of (18) in the following form 

a (z, t) = yi (t) exp (iQz) , b (z, t) = y 2 (t) exp (iQ.z) , 

where y\ (t) et y 2 (t) are two functions and is an arbitrary constant. Then 
we obtain the system 

m + {vi + yl) yi = ( fi - ^o) yi, 
m + (yl + yl) y 2 = {^i + %) yi- 

The latter coincides obviously with (16) for Ai = Q — fio and A2 = Q + fio- 
2.6 The Yang-Mills equations. 

We consider the Yang-Mills system for a field with gauge group SU (2) : 

where F jk ,Aj G T e SU(2),l < j,k < 4 and F jk = ^ - ^ + [A,-,A fc ]. 
The self-dual Yang-Mills (SDYM) equations is an universal system for which 
some reductions include all classical tops from Euler to Kowalewski (0+1- 
dimensions), K-dV, Nonlinear Schrodinger, Sine-Gordon, Toda lattice and N- 
waves equations (1+1-dimensions), KP and D-S equations (2+1-dimensions). 
In the case of homogeneous double-component field, we have djA k = 0, j 7^ 
l,Ai = A 2 = 0, A 3 = mUi G su(2),A 4 = n 2 U 2 G su(2) where n { are s«(2)- 
generators (i.e., they satisfy commutation relations : n\ = [n 2 , [ni,n 2 ]],n 2 = 
[ni, [n 2 ,ni]}). The system becomes 

-^ + ^2=0, -^ + Wl 2 = 0, 

with t = r\. By setting ?7j = (/j, = pj, j = 1, 2, Yang-Mills equations are 
reduced to hamiltonian system 

• 7^ r u 7 f - 1 

X= ~dx' X = ^ 1,92 ' Pl ' P2 ) ' 1 J O 

with i7 = ^(Pi+p 2 + q 2 q2), the hamiltonian. The symplectic transformation 
Pi ^ ^{pi+P2),P2 -r\ ^-(p\-p<i),q\ -r\ ±(\/2)(q 1 +iq 2 ),q 2 ^ 5(^2)^1- 
152); takes this hamiltonian into 

# = ^(p? + pi) + + + l<ihl 

which coincides with (15) for Ai = A2 = 0. 
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3 Algebraic complete integrability 



We give some results about abelian surfaces which will be used, as well as the 
basic techniques to study two-dimensional algebraic completely integrable 
systems. Let M = C/A be a n— dimensional abelian variety where A is the 
lattice generated by the In columns Ai, . . . , \2n of the n x In period matrix 
O and let D be a divisor on M. Define C(T>) = {/meromorphic onM : (/) > 
—V}, i.e., for V = ^ kjVj a function / G CiV) has at worst a kj— fold pole 
along Vj. The divisor V is called ample when a basis (/o, . . . , /at) of C(kV) 
embeds M smoothly into P^ for some k, via the map M — > P^, p w [1 : 
flip) ■ ... : /jv(f)]> then kV is called very ample. It is known that every 
positive divisor T> on an irreducible abelian variety is ample and thus some 
multiple of V embeds M into P^. By a theorem of Lefschetz, any k > 3 
will work. Moreover, there exists a complex basis of C n such that the lattice 
expressed in that basis is generated by the columns of the n x 2n period 
matrix 

/ <*i | \ 

| z . 

V o s n | / 

with Z T = Z,lmZ > 0, 6j G N* and 5j|(5j + i. The integers <5j which provide 
the so-called polarization of the abelian variety M are then related to the 
divisor as follows : 

dim£(D) = Si ...S n . (19) 

In the case of a 2— dimensional abelian varieties (surfaces), even more can be 
stated : the geometric genus g of a positive divisor V (containing possibly 
one or several curves) on a surface M is given by the adjunction formula 

K M -V + V.V , 1 fnn . 

gi v ) = g + ' ( 2 °) 

where Km is the canonical divisor on M, i.e., the zero-locus of a holomorphic 
2— form, T>.T> denote the number of intersection points of T> with a+T> (where 
a + V is a small translation by a of V on M), where as the Riemann-Roch 
theorem for line bundles on a surface tells you that 

xiV)= Pa iM) + l + ±iV.V-VK M ), (21) 

where p a (M) is the arithmetic genus of M and x(P) the Euler characteristic 
of V. To study abelian surfaces using Riemann surfaces on these surfaces, 
we recall that 

X(D) = dimH iM,O M iV)) - dim F X (M, M p)), 

= dim£(£>) - dim H^M, ft 2 (£> K* M )), (Kodaira-Serre duality), 
= dim CiV), (Kodaira vanishing theorem), (22) 
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whenever V <g) K* M defines a positive line bundle. However for abelian sur- 
faces, Km is trivial and p a (M) = — 1; therefore combining relations (19), 
(20), (21) and (22), 

X (V) = dim£(P) = ^= g (V)-l = 5^2. 

A divisor T> is called projectively normal, when the natural map C{V)® k — ► 
C(kV), is surjective, i.e., every function of C(kV) can be written as a linear 
combination of k-fold products of functions of C(T>). Not every very ample 
divisor T> is projectively normal but if T> is linearly equivalent to kT>o for 
k > 3 for some divisor Vq, then V is projectively normal. 
Now consider the exact sheaf sequence 

— ► O c O g — ► ^ — ► 0, 

where C is a singular connected Riemann surface, C = Cj the correspond- 
ing set of smooth Riemann surfaces after desingularization and ir : C — > C 
the projection. The exactness of the sheaf sequence shows that the Euler 
characteristic 

X(P) = dimH (O) - dimH^O), 

satisfy 

X(O c ) ~ X{0~) + X(X) = 0, (23) 

where X(X) only accounts for the singular points p of C; X(X p ) is the di- 
mension of the set of holomorphic functions on the different branches around 
p taken separately, modulo the holomorphic functions on the Riemann sur- 
face C near that singular point. Consider the case of a planar singularity (in 
this paper, we will be concerned by a tacnode for which X(X) = 2, as well), 
i.e., the tangents to the branches lie in a plane. If fj(x,y) = denote the 
j th branch of C running through p with local parameter Sj, then 

*,x p ) = di m n 3 c[ N ]/||i. 

So using (22) and Serre duality, we obtain X(Oc) = 1 — g(C) and X(Og) = 
n — ^j=i5 f (Ci)- Also, replacing in the formula (23), gives 

n 

g{C)=Y J 9(C 3 ) + X{X) + l-n. 

3=1 

Finally, recall that a Kahler variety is a variety with a Kahler metric, i.e., a 
hermitian metric whose associated differential 2-form of type (1, 1) is closed. 
The complex torus C 2 /lattice with the euclidean metric ^2dzi <g> dzi is a 
Kahler variety and any compact complex variety that can be embedded in 
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projective space is also a Kahler variety. Now, a compact complex Kahler 
variety having as many independent meromorphic functions as its dimension 
is a projective variety 

Consider now hamiltonian problems of the form 

BH 

X H :x = J— = f(x), x G K m , (24) 

where H is the hamiltonian and J = J(x) is a skew-symmetric matrix 
with polynomial entries in x, for which the corresponding Poisson bracket 
{Hi,Hj} = (^N^^r)) satisfies the Jacobi identities. The system (24) 
with polynomial right hand side will be called algebraic complete integrable 
(a.c.i.) when : 

i) The system possesses n+k independent polynomial invariants H±, H n+ ^ 

OH 

of which k lead to zero vector fields J d n x +l (x) = 0, 1 < i < k, the n re- 
maining ones are in involution (i.e., {Hi, Hj} = 0) and m = 2n + /c. For most 

n+k 

values of Cj G M, the invariant varieties f] {x £ M. m : Hi = Cj} are assumed 

i=i 

compact and connected. Then, according to the Arnold-Liouville theorem, 
there exists a diffeomorphism 

n+k 

p| {x e M m : Hi = Ci} -► R n / Lattice, 
i=l 

and the solutions of the system (24) are straight lines motions on these tori. 

ii) The invariant varieties, thought of as affine varieties in C m can be com- 
pleted into complex algebraic tori, i.e., 

n+k 

p| {Hi = Ci,xe C m } UV = C n /Lattice, 
i=i 

where C n / Lattice is a complex algebraic torus (i.e., abelian variety) and V 
a divisor. Algebraic means that the torus can be defined as an intersection 

M 

Pi(Xo, Xn) = 0} involving a large number of homogeneous polynomi- 
al 

als Pj. In the natural coordinates (ti, ...,t n ) of C n /Lattice coming from C n , 
the functions X{ = Xi{t\, t n ) are meromorphic and (24) defines straight 
line motion on C n / Lattice. Condition i) means, in particular, there is an al- 
gebraic map {x\(t), ...,x m (t)) i— > (fii(t), ...,[x n (t)) making the following sums 
linear in t : 




= djt , 1 < j < n, dj G C, 



where oj\,...,u n denote holomorphic differentials on some algebraic curves. 
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The existence of a coherent set of Laurent solutions 



oo 

Xi = j^a^V - **, ki G Z, some h > 0, (25) 
3=0 

depending on dim (phase space) — 1 = m — 1 free parameters is necessary 
and sufficient for a hamiltonian system with the right number of constants 
of motion to be a.c.i. So, if the hamiltonian flow (24) is a.c.i., it means 
that the variables Xj are meromorphic on the torus C n / Lattice and by com- 
pactness they must blow up along a codimension one subvariety (a divisor) 
V C C™ / 'Lattice. By the a.c.i. definition, the flow (24) is a straight line 
motion in C n / Lattice and thus it must hit the divisor V in at least one 
place. Moreover through every point of T>, there is a straight line motion 
and therefore a Laurent expansion around that point of intersection. Hence 
the differential equations must admit Laurent expansions which depend on 
the n — 1 parameters defining V and the n + k constants Cj defining the torus 
C n /Lattice , the total count is therefore m — 1 = dim (phase space) — 1 
parameters. 

Assume now hamiltonian flows to be (weight)-homogeneous with a weight 
Vi G N, going with each variable x«, i.e., 

fi (a ui xi, ...,a Um x m ) = a Ut+1 fi (x±, ...,x m ) , Vq G C. 

Observe that then the constants of the motion H can be chosen to be 
(weight)-homogeneous : 

H (a Ul xi, ...,a Um x m ) = a k H (x±, ...,x m ) , k G Z. 

If the flow is algebraically completely integrable, the differential equations 
(24) must admits Laurent series solutions (25) depending on m — 1 free 
parameters. We must have ki = Vi and coefficients in the series must satisfy 
at the th step non-linear equations, 

fi (xf\ ..,x&)) + 9l xf ] = 0, 1 < i < m, (26) 
and at the k* ft step, linear systems of equations : 

{L-kI)zW = l , . , ° % k = \ k 1)f u , (27) 

[ some polynomial m x v x y ' tor k > 1, 

where 

L = Jacobian map of (26) = — — h gl L_.(o) • 

oz 

If m — 1 free parameters are to appear in the Laurent series, they must either 
come from the non-linear equations (26) or from the eigenvalue problem (27) , 
i.e., L must have at least m — 1 integer eigenvalues. These are much less 
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conditions than expected, because of the fact that the homogeneity k of the 
constant H must be an eigenvalue of L Moreover the formal series solutions 
are convergent as a consequence of the majorant method. Next we assume 
that the divisor is very ample and in addition projectively normal. Consider 
a point p £ V, a chart Uj around p on the torus and a function yj in C(V) 
having a pole of maximal order at p. Then the vector (1/yj, yi/yj, ■ ■ ■ , UN/Vj) 
provides a good system of coordinates in Uj. Then taking the derivative with 
regard to one of the flows 

,My = ViVj-ViVj l<j<N, 

are finite on Uj as well. Therefore, since y| has a double pole along V, the 
numerator must also have a double pole (at worst), i.e., yiyj — y%yj £ C(2V). 
Hence, when V is projectively normal, we have that 

(— =i^ a k,i{— )(— )> 

Vi Xi Vj Vj 

i.e., the ratios yi/yj form a closed system of coordinates under differentiation. 
At the bad points, the concept of projective normality play an important role: 
this enables one to show that yi/yj is a bona fide Taylor series starting from 
every point in a neighbourhood of the point in question. 

To prove the algebraic complete integrability of a given hamiltonian sys- 
tem, the main steps of the method are : 

- The first step is to show the existence of the Laurent solutions, which 
requires an argument precisely every time k is an integer eigenvalue of L and 
therefore L — kl is not invertible. 

- One shows the existence of the remaining constants of the motion in 
involution so as to reach the number n + k. 

- For given Ci, ...,c m , the set 



V = 



(0) , _ L „(2) +2 

i 

Laurent solutions such that : Hj (xi (t)) = Cj + Taylor part 



Xi (t) = t~ Ui ( x\ > + x\ 't + x\ V + ■ ■ • J , 1 < i < m, 

Ij (Xi (t)) = Cj 



defines one or several n — 1 dimentional algebraic varieties (divisor) having 
the property that 

n+k 

P| {Hi = a,z G C m } UV = a smooth compact, connected variety 



=i 



with n commuting vector fields 
independent at every point. 
= a complex algebraic torus T n = C n / Lattice. 
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The flows J 



dH, 



k + i 



J 



dH, 



k-\-n 



are straight line motions on T n . 



dz ' "■' dz 

From the divisor V, a lot of information can be obtained with regard to the 
periods and the action-angle variables. 



3.1 A five- dimensional system 

Consider the following system of five differential equations in the unknowns 
z 1 ,...,z 5 : 



zi = 2z 4 , z 3 = z 2 {3z! + 8z 2 ), 

z 2 = z 3 , i 4 = z\ + kz\z\ + z 5 , 



(28) 



Z5 = 2z\ Z4 + 4^2-^4 — 2z\z 2 z 3 . 
The following three quartics are constants of motion for this system 

F i = ~ ZlZ % + 2 Z * ~ 1 Z * ~ 2Z *' 

F 2 = z^ - zfz 5 + 4z 1 z 2 z 3 z 4 - zfzj + ^zf - 4z 2 zf, 



(29) 



^3 



■ 2 2 
ZiZ 5 + z 1 z 2 



This system is completely integrable and the hamiltonian structure is de- 
fined by the Poisson bracket {F, H} = (§f , J™) = ^^ i=1 4*gf§, where 

dH _ { dH dH_ dH_ dH_ dH_\T QT1 j 

^zi ' dz 2 ' 9z 3 ' d Zi > dz 5 > ' dI1Q 



J 





1 
0-10 

-2zi 

-4z 4 4zi z 2 



2zi 




-2z 5 + 8^12:2 



4z 4 


~4:ZlZ 2 

2z§ — %z\z\ 




is a skew-symmetric matrix for which the corresponding Poisson bracket sat- 
isfies the Jacobi identities. The system (28) can be written as z = J^jr, z = 
(z±, z 2 , z 3 , Z4, z^) T , where H = F\. The second flow commuting with the first 
is regulated by the equations z = J^r, z = (z±, z 2 , z 3 , 24, z^) T . These vector 
fields are in involution : {F±, F 2 } = (^ L , J^§r) = 0> an d the remaining one 



is casimir : J 



dz 



0. The invariant variety A defined by 



A = f]{z : F k (z) = c k } C 



(30) 



k=i 



is a smooth affine surface for generic values of (ci, c 2 , c 3 ) G C 3 . So, the ques- 
tion I address is how does one find the compactification of A into an abelian 
surface? The idea of the direct proof we shall give here is closely related 
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to the geometric spirit of the (real) Arnold- Liouville theorem. Namely, a 
compact complex n- dimensional variety on which there exist n holomorphic 
commuting vector fields which are independent at every point is analytically 
isomorphic to a n-dimensional complex torus C n / Lattice and the complex 
flows generated by the vector fields are straight lines on this complex torus. 
Now, the main problem will be to complete A (30) into a non singular com- 
pact complex algebraic variety A = A U V in such a way that the vector 
fields Xp 1 and Xp 2 generated respectively by F\ and F2, extend holomor- 
phically along a divisor V and remain independent there. If this is possible, 
A is an algebraic complex torus (an abelian variety) and the coordinates 
zi, . . . , £5 restricted to A are abelian functions. A naive guess would be to 
take the natural compactification A of A by projectivizing the equations: 
A = ff k =i{ F k( z ) = Cfc^ol c p5 - Indeed, this can never work for a general 
reason : an abelian variety A of dimension bigger or equal than two is never 
a complete intersection, that is it can never be described in some projective 
space P n by n-dim A global polynomial homogeneous equations. In other 
words, if A is to be the affine part of an abelian surface, A must have a 
singularity somewhere along the locus at infinity A n {Zq = 0} . In fact, we 
shall show that the existence of meromorphic solutions to the differential 
equations (28) depending on 4 free parameters can be used to manufacture 
the tori, without ever going through the delicate procedure of blowing up 
and down. Information about the tori can then be gathered from the divisor. 

Theorem 10 The system (28) possesses Laurent series solutions which de- 
pend on 4 free parameters : a, j3, 7 and 0. These meromorphic solutions 
restricted to the surface A(30) are parameterized by two copies C_i and C\ 
of the same Riemann surface (32) of genus 7. 

Proof. The first fact to observe is that if the system is to have Laurent solu- 
tions depending on 4 free parameters a, (3, 7, 9, the Laurent decomposition 
of such asymptotic solutions must have the following form 

zi = -a- \o? +pt-—a (a 3 + 4/?) t 2 + 7 t 3 + • • • , 

z 2 = Y t £ ~ \ ea + l £a2f ~ (-" 3 + 12 ^) t 2 + et 3 + -- - , 

z 3 = -^e + ^ea 2 -^e(-a 3 + 12p)t + 39t 2 + --- , (31) 
za = ~^2 a + \P- ^a(a 3 + 4/3)t + ^ 7 t 2 + --- , 

z 5 = ^ a 2 - (a 3 + 4/3) + (a 3 + 2(3) - (a 2 (3 - 2 7 + 4e(9a) * + ■■■, 

with e = ±1. Using the majorant method, we can show that these series 
are convergent. Substituting the Laurent solutions (31) into (29): F\ = c\, 
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F2 = C2 and F% = c%, and equating the t°-terms yields 

7 15 

Fi = —a 4 - -a(3 - -eO = a, 

F 2 = ^ (4/3-a 3 ) (4a 2 /3-a 5 + 64^a-32 7 ) =c 2 , 
F 3 = --^a 6 - (3 2 -ja 3 (3 -3e6a 2 + 4a~f = c 3 . 

Eliminating 7 and 6 from these equations, leads to an equation connecting 
the two remaining parameters a and (3 : 

C : 64/3 3 - 16a 3 /? 2 - 4 (a 6 - 32a 2 ci - 16c 3 ) 

+a (32c 2 - 32q 4 ci + a 8 - 16a 2 c 3 ) = 0. (32) 

The Laurent solutions restricted to the surface A(30) are thus parameter- 
ized by two copies C-\ and C\ of the same Riemann surface C(32). According 
to the Riemann-Hurwitz formula, the genus of the Riemann surface C is 7, 
which establishes the theorem. 



In order to embed C into some projective space, one of the key underly- 
ing principles used is the Kodaira embedding theorem, which states that a 
smooth complex manifold can be smoothly embedded into projective space 
P N with the set of functions having a pole of order k along positive di- 
visor on the manifold, provided k is large enough; fortunately, for abelian 
varieties, k need not be larger than three according to Lefshetz. These func- 
tions are easily constructed from the Laurent solutions (31) by looking for 
polynomials in the phase variables which in the expansions have at most a 
k-fold pole. The nature of the expansions and some algebraic proprieties of 
abelian varieties provide a recipe for when to terminate our search for such 
functions, thus making the procedure implementable. Precisely, we wish to 
find a set of polynomial functions {/o, . . . , /at}, of increasing degree in the 
original variables z\,...,z^ having the property that the embedding T> of 
C\ +C_i into P N via those functions satisfies the relation : geometric genus 
(V) = g(T>) = N + 2. A this point, it may be not so clear why V must really 
live on an abelian surface. Let us say, for the moment, that the equations of 
the divisor V (i.e., the place where the solutions blow up), as a Riemann sur- 
face traced on the abelian surface A (to be constructed in theorem 12), must 
be understood as relations connecting the free parameters as they appear 
firstly in the expansions (31). This means that (32) must be understood as 
relations connecting a and j3. Let 



L (r) 



polynomials / = f(z t . . . , 25) 
of degree < r, such that 
f{z{t))=t-\z^ +...), 
with / on V 

and with z(t) as in (4) 



\/[F k = c k ,k = 1,2,3], 
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and let (fa, fa, ... , /iv r ) be a basis of L^. We look for r such that : g(T>^) = 
N r + 2,2?( r ) C ¥ Nr . We shall show that it is unnecessary to go beyond r=4. 



Theorem 11 a) The spaces L^', nested according to weighted degree, are 
generated as follows 



L« = 


{fa, fa, fa}, 




= 


L {1) @{fa,h,fa,fa], 




= 


L {2) @{h,fa,fa,f W }, 




LW = 


^ (3) ©Ul2,/l3,/l4,/l 5 }, 


(33) 


zi,fa = 


Z2J3 = 2z 5 - z\,fa = z 3 + 2ez\,fa 


= Z4 + 



ez lZ2 J 6 = [/i,/ 2 ] ,/ 7 = + 2e/ 4 ),/ 8 = fa(fi + 2e/ 4 ),/ 9 = z 4 (/ 3 + 

2efa), /10 = 2 5 (/ 3 + 2e/ 6 ), /11 = / 5 (/i + 2e/ 4 ), /12 = /i/ 2 (/ 3 + 2efa), /13 = 
/4/s + A] , /i4 = / 3 ] + 2e [fa, fa] , fi 5 = fa ~ 2z 5 + 4/|, unto fo, a fc ] = 
SjSfc — SjSfc, toe wronskien of Sk and Sj. 

b) provides an embedding ofD^ into projective space P 15 andV>^ has 
genus 17. 



Proof, a) The proof of a) is straightforward and can be done by inspection 
of the expansions (31). 

b) It turns out that neither nor L(2) , nor L^ 3 ) , yields a Riemann surface 
of the right genus; in fact g(T>^) 7^ dimL^ + l,r = 1,2,3. For instance, 
the embedding into P 2 via does not separate the sheets, so we pro- 
ceed to and the corresponding embedding into P 6 is unacceptable since 
5 (£>(2)) - 2 > 6 and C P 6 + P 9 " 2 , which contradicts the fact that 
N r = g(V^) — 2. So we proceed to and we consider the correspond- 
ing embedding into P 10 , according to the functions (fo, . . . , fao). For finite 
values of a and fi, dividing the vector (/o,---,/io) by fa and taking the 
limit t -> 0, to yield [0 : 2ea : 1 : -e(Afi - a 3 ) : -a : -ea 2 : §(4/3 - a 3 ) : 
ea 3 : §a 2 : |ea 3 (4/3 — a 3 ) : — \eofi(Afi — a 3 )]. The point a = require 
special attention. Indeed near a = 0, the parameter /3 behaves as follows 
: fi ~ 0, iy/c~3~, — iy/c~3~. Thus near (a,fi) = (0,0), the corresponding point 
is mapped into the point [0 : : 1 : : : : : : : : 0] in P 10 
which is independent of e = ±1, whereas near the point (a,fi) = (0, i-s/c-s) 
(resp. (ct,fi) = (0,—iy/c3)) leads to two different points : [0:0:1: 
-Aei^fci : : : 2eijc^ : : : : 0] (resp. [0:0:1: Aei^i : : 
: — 2eiy/c~3 : : : : 0]), according to the sign of e. The Riemann sur- 
face (31) has three points covering a = 00, at which fi behaves as follows 

: P W° 3 ' ik" 3 ( 1333 ~ 1295i>/3) , ^a 3 (1333 + 1295i^/3~) . Then by 

dividing the vector (fa, . . . , fao) by fao, the corresponding point is mapped 
into the point [0:0:0:0:0:0:0:0:0:0:1] in P 10 . Thus, 
g (vW) - 2 > 10 and C P 10 / P s ~ 2 , which contradicts the fact that 
N r = g(V^)) - 2. Consider now the embedding into P 15 using the 16 
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functions /o, . . . , /15 of L^(33). It is easily seen that these functions separate 
all points of the Riemann surface (except perhaps for the points at a = 00 
and a = f3 = 0) : The Riemann surfaces C\ and C-\ are disjoint for finite 
values of a and (3 except for a = (3 = 0; dividing the vector (/o, . . . , /15) 
by /2 and taking the limit t — > 0, to yield [0 : 2ea : 1 : —e(4(3 — a 3 ) : 
-a : -ea 2 : ±(4/3 - a 3 ) : ea 3 : \c? : \ea 3 {4(3 - a 3 ) : -±ea 4 (4/3 - a 3 ) : 
-±ea 4 : -±a 3 (4/3 - a 3 ) : fa (4/3 - a 3 ) : ea 3 (4/3 - a 3 ) : -2ea 3 }, As before, 
the point a = require special attention and the parameter (3 behaves as 
follows : (3 ~ 0, iy/cH, — iy/cii. Thus near (a, (3) = (0,0), the corresponding 
point is mapped into the point [0:0:1:0:0:0:0:0:0:0:0:0: 
: : : 0] in P 15 which is independent of e = ±1, whereas near the point 
(a, /?) = (0, iy/ci) (resp. (a, (3) = (0, — iy/cH)) leads to two different points 
: [0:0:1: -Aei^ : : : lei^fci :0:0:0:0:0:0:0:0:0] (resp. 
[0:0:1: 4ei^ : : : -2ei^ci :0:0:0:0:0:0:0:0: 0]), 
according to the sign of e. About the point a = 00, it is appropriate 
to divide by /10; then the corresponding point is mapped into the point 
[0 : : : : : : : : : : 1 : : : : : 0], in P 15 which is 
independent of e. The divisor T>^ obtained in this way has genus 17 and 
p( 4 ) p 15 = ps~ 2 ; as desired. This ends the proof of the theorem. 

Let C = L (4) and V = £> (4) . Next we wish to construct a surface strip 
around V which will support the commuting vector fields. In fact, V has 
a good chance to be very ample divisor on an abelian surface, still to be 
constructed. 

Theorem 12 The variety A (30) generically is the affine part of an abelian 
surface A. The reduced divisor at infinity A\A = C\ + C_i, consists of two 
copies C\ and C-i of the same genus 7 Riemann surface C(32). The system 
of differential equations (28) is algebraically completely integrable and the 
corresponding flows evolve on A. 

Proof. We need to attaches the affine part of the intersection of the three 
invariants ^1,^2,^3 so as to obtain a smooth compact connected surface 
in P 15 . To be precise, the orbits of the vector field (28) running through 
T> form a smooth surface S near T> such that T,\A C A and the variety 
A = Au X is smooth, compact and connected. Indeed, let tp(t,p) = {z(t) = 
(zi(t), z 5 (t)) : t £ C, < \t\ < e}, be the orbit of the vector field (28) go- 
ing through the point p £ A. Let S p C P 15 be the surface element formed by 
the divisor V and the orbits going through p, and set E = U pe x>E p . Consider 
the Riemann surface V = TL PI E where TL C P 15 is a hyperplane transversal 
to the direction of the flow. If V is smooth, then using the implicit function 
theorem the surface E is smooth. But if V is singular at 0, then E would be 
singular along the trajectory (t— axis) which go immediately into the affine 
part A. Hence, A would be singular which is a contradiction because A is 
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the fibre of a morphism from C 5 to C 3 and so smooth for almost all the 
three constants of the motion c&. Next, let A be the projective closure of 
A into P 5 , let Z = [Z : Z x : . . . : Z 5 ] G P^_and let I = A n {Z = 0} 
be the locus at infinity. Consider the map A C P 5 — > P 15 , Z t— > f(Z), 
where / = (/o, /i, /15) G £(1?) and let ^4 = /(A). In a neighbourhood 
V(p) C P 15 of p, we have S p = j4 and C A. Otherwise there would 

exist an element of surface T,' p C A such that S p n Ej, = (t - axis), or- 
bit tp(t,p) = (t — axis)\ p C A, and hence A would be singular along the 
t— axis which is impossible. Since the variety A n {Zo 7^ 0} is irreducible 
and since the generic hyperplane section Ti. gen . of A is also irreducible, all 
hyperplane sections are connected and hence I is also connected. Now, con- 
sider the graph T f C P 5 x P 15 of the map /, which is irreducible together 
with A. It follows from the irreducibility of I that a generic hyperplane 
section Tf n {Ti. gen . x P 15 } is irreducible, hence the special hyperplane sec- 
tion Tf n {{Zq = 0} x P 15 } is connected and therefore the projection map 
proj ¥ ir,{Tf n {{Z Q = 0} x P 15 }} = /(/) = V, is connected. Hence, the 
variety A U £ = A is compact, connected and embeds smoothly into P 15 
via /. We wish to show that A is an abelian surface equipped with two ev- 
erywhere independent commuting vector fields. For doing that, let <p T1 and 
4> T2 be the flows corresponding to vector fields Xp 1 and Xp 2 . The latter 
are generated respectively by Fi and F2. For p G V and for small e > 0, 
4> T1 (p) , Vri , < I Ti| < e, is well defined and (j) T1 (p) G A\A. Then we may de- 
fine cP on A by 4> T2 (q) = (f)- Tl (f) T2 (f) Tl (q),q G U(p) = </r n (U((f) Tl (p))), where 
f7 (p) is a neighbourhood of p. By commutativity one can see that 4> T2 is in- 
dependent of n; c/)- Tl - £l (t) T2 cl) Tl+£l (q) = cf)- Tl (t)- £l (t) T2 (j) T1 (t) £l = (p- Tl 4> T2 (p Tl (q). 
We affirm that 4> T2 {q) is holomorphic away from V. This because (p T2 (j) T1 (q) 
is holomorphic away from V and that (j) Tl is holomorphic in U (p) and maps 
bi-holomorphically U(p) onto U(cp T1 (p)). Now, since the flows (/> ri and 4> T ' 2 are 
holomorphic and independent on V, we can show along the same lines as in 
the Arnold-Liouville theorem [15] that A is a complex torus C 2 / lattice and 
so in particular A is a Kahler variety. And that will done, by considering the 
local diffeomorphism C 2 — > A, (ti,T2) 1— > <fi Tl <ft T2 (p) , for a fixed origin p £ A. 
The additive subgroup {(ti,T2) G C 2 : (p T1 (p T2 (p) = p} is a lattice of C 2 , 
hence C 2 /lattice —> A is a biholomorphic diffeomorphism and ^4 is a Kahler 
variety with Kahler metric given by dr\ <g) cffi + c?T2 <8> c?T2- As mentioned in 
appendix A, a compact complex Kahler variety having the required number 
as (its dimension) of independent meromorphic functions is a projective va- 
riety. In fact, here we have A C P 15 . Thus A is both a projective variety and 
a complex torus C 2 /lattice and hence an abelian surface as a consequence 
of Chow theorem. This completes the proof of the theorem. 

Remark 3.1 a) Note that the reflection a on the affine variety A amounts to 
the flip a : (zi, Z2, £3, 24, £5) ► (zi, — Z2, £3, — £4, 25), changing the direction 
of the commuting vector fields. It can be extended to the (-Id) -involution 
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about the origin of C 2 to the time flip (t±,t2) > (— ii,— £2) on A, where 
ti and ti are the time coordinates of each of the flows Xp l and Xp 2 . The 
involution a acts on the parameters of the Laurent solution (30) as follows a : 
(t,a,/3,j,9) 1 — ► (— t, — a, — /3, — 7, 9), interchanges the Riemann surfaces C £ 
and the linear space C can be split into a direct sum of even and odd functions. 
Geometrically, this involution interchanges C\ and C_i, i.e., C_i = oC\. 
b) Consider on A the holomorphic 1-forms dt\ and a%2 defined by dti(Xp j ) = 
5ij, where Xp ± and Xp 2 are the vector fields generated respectively by F\ and 
F2. Taking the differentials of C, = l/z\ and £ = Z1/Z2 viewed as functions 
of t\ and t2, using the vector fields and the Laurent series (31) and solving 
linearly for dt\ and dt2, we obtain the holomorphic differentials 

1 ,d£ ,, d( 8 
Ul = dh ^ = A ( ^ C "^ )k = a (-40 -Ho*)**' 



^ = d ^ = ~K ~dii ~ ~dt[ = (-40 + 03)'**' 

with A = — JlrJ^"- The zeroes of u>2 provide the points of tangency of 

the vector field Xp 1 to C £ . We have ^ = — (—4/3 + a 3 ) , and Xp 1 is tangent 
to 7i £ o,t the point covering a = 00. 

3.2 The Henon-Heiles system 

The Henon-Heiles system 

dH OH 8H dH 

qi = -k— , Q2 = -5—, pi = , P2 = , (34) 

dpi dp 2 <% 

with 

H = H 1 = -(p 2 +p 2 + aq 2 + 6g 2 ) + g 2 o 2 + 6^, 

has another constant of motion 

H 2 = q\ + ^q\q\ - 4pi (pig 2 - f>29i) + 4ag 2 a 2 + (4a - b) (pf + aqf) , 

where a, 6, are constant parameters and qi, q2,Pi,P2 are canonical coordi- 
nates and momenta, respectively. First studied as a mathematical model 
to describe the chaotic motion of a test star in an axisymmetric galactic 
mean gravitational field this system is widely explored in other branches of 
physics. It well-known from applications in stellar dynamics, statistical me- 
chanics and quantum mechanics. It provides a model for the oscillations of 
atoms in a three-atomic molecule. The system (34) possesses Laurent series 
solutions depending on 3 free parameters a, (3, 7, namely 

a fa 3 aA aB\ r2 , U) ,3 , J5) ,4 , (6). 5 , 
1 | a 2 B t fa 4 a 2 A a 2 B B 2 \^ 2 i a/3 3 4 



q2 = -¥ + T2-n + \l8 + ^--m-^6) t+ T t+lt + 
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where p\ = q l7 p 2 = k 2 an d 



,( 4 ) 



aAB a 5 lla 3 B lla 3 A aB 2 aA 2 



qi 24 72 + 720 120 720 8 ' 

(5 ) _pa^_ f3B_A£ 
Ql 12 60 10 ' 

(6 ) a 7 a 1 a 5 A a b B a 3 B 2 a 3 A 2 a 3 AB aAB 2 
Ql 9~ ~ 15552 ~ 2160 + 12960 + 25920 + 1440 ~ 4320 + 1440 

aB 3 aA 2 B aA 3 
~ 19440 ~ 288 + "144 ' 

Let V be the pole solutions restricted to the surface 

2 

M c = P| {x = {qi,q2,Pi,P2) e C 4 ,Hi (x) = q} , 
i=i 

to be precise V is the closure of the continuous components of the set of 
Laurent series solutions x (t) such that iJj(x(i)) = Cj, 1 < i < 2, i.e., 
V = t° — coefficient of M c . Thus we find an algebraic curve defined by 

V : (3 2 = P 8 (a), (35) 

where 

8V ; 15552 432 V 18 / 36 V 15120 7 1260 ) 

1 / 1 , 13 j9 2 , 9 10 \ o 1 

4A 3 B 3 A 2 B + -AB 2 - — Cl a 2 + — c 2 . 

36 V 2520 6 9 7 J 36 

The curve V determined by an eight-order equation is smooth, hyperelliptic 
and its genus is 3. Moreover, the map 

a:V^V, (/?,<*)■— (/?,-«), (36) 

is an involution on V and the quotient £ = V/a is an elliptic curve defined 
by 

£ : P 2 = P 4 (C), (37) 

where P4 (£) is the degree 4 polynomial in C = ot 2 obtained from (35) . The 
hyperelliptic curve V is thus a 2-sheeted ramified covering of the elliptic 
curve £ (37) , 

p:V^£, (fta).— (/J,C), (38) 
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ramified at the four points covering £ = and 00. The affine surface M ( 
completes into an abelian surface M c , by adjoining the divisor V. The latter 
defines on M c a polarization (1,2). The divisor 2T> is very ample and the 
functions 1, yi, yf, y 2 , 21, xj + y 2 y 2 , x 2 yi~2x 1 y 2 , x 1 x 2 + 2Ay 1 y 2 + 2y 1 y%, 
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embed M c smoothly into CP 7 with polarization (2,4). Then the system (34) 
is algebraically completely integrable and the corresponding flow evolues on 
an abelian surface M c = C 2 /lattice, where the lattice is generated by the 

period matrix ( 2 Q \ « \ ) , Im (° ^ > 0. 

Theorem 13 The abelian surface M c which completes the affine surface M c 
is the dual Prym variety Prym* (T>/£) of the genus 3 hyperelliptic curve T> 
(35) for the involution a interchanging the sheets of the double covering p 
(38) and the problem linearizes on this variety. 

Proof. Let (01,02,03,61,62,63) be a canonical homology basis of V such 
that cr(cti) = 03, cr(&i) = 63, a (02) = —02, a (62) = —62, for the in- 
volution a (36). As a basis of holomorphic differentials luq,uji,uj2 on the 
curve V (35) we take the differentials wi = ^,w 2 = 7^3 = 2s*s, 
and obviously a*(uj\) = — u\, a*(u2) = —U2,-(t*(uj^) = u 3 . Recall that the 
Prym variety Prym(V/£) is a subabelian variety of the Jacobi variety 
Jac{V) = Pic°(V) = H l {O v ) / H l (V,Z) constructed from the double 
cover p : the involution a on T> interchanging sheets, extends by linearity 
to a map a : JaciV) — > JaciV) and up to some points of order two, Jac(V) 
splits into an even part and an odd part : the even part is an elliptic curve 
(the quotient of V by a, i.e., £ (18)) and the odd part is a 2— dimensional 
abelian surface Prym (V /£) . We consider the period matrix Q of Jac(D) 



Then, 



/ 4 Wl Ia 2 UJl fa 3 U l h U l Ib 2 Ul fbs U l 
= 4^ J a2 "2 4^2 f bl U2 4^2 j b ^2 
\ la, ^3 J a2 ^3 J a3 UJ 3 J bi UJ 3 J b2 UJ 3 ^ W 3 



= I / ai ^2 4 ^2 - 4 UJ 2 4 W 2 4 W 2 - 4 W 2 
4 W3 iai^ 3 Al^ 3 ^i^ 3 

and therefore the period matrices of Jac(£) (i.e., £), Prym(V / £) and Prym*(V/£) 
are respectively A = (4 ^3 4 ^ 3 )' 

2 4^ 2 24^2 4^2 



and 



4 Wi 4 Wi 4 Wi 
4^2 4^2 4^2 4^2 
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u 2 +n.i f b . [ uj 2 J : mi,m € Z}, be the period 

U 3 J \ U 3 

lattice associated to fi. Let us denote also by La, the period lattice associated 
A. We have the following diagram 



I 

£ V 

IV* / if 

— ► kerA^ — ► Prym(V /£ © £ = Jac{V) ^ £ — > 

\r _ | 

M c = M c U 2V ~ C 2 /lattice 

I 


The polarization map r : Prym(V/£) — ► M c = Prym* (V/£), has kernel 
(<f*£ ) ~ Z2XZ2 and the induced polarization on PrymiV j '£) is of type (1,2). 
Let M c — ► C 2 /L\ : p rx (dtl)' ^ e ^ ne uniformizing map where dii, ^2 are 
two differentials on M c corresponding to the flows generated respectively by 
H\,H 2 such that : dti\x> = w\ and dt 2 \v = u> 2 , 



fc=l \ JV k 

is the lattice associated to the period matrix 

L,dti f U2 dh f U4 dh f U4 dh 
J Ul dt 2 J u2 dt 2 J u3 dt 2 J u4 dt 2 

and (i>\, z^2, ^3, fi) is a basis of Hi(M c , Z). By the Lefschetz theorem on hy- 
perplane section [9], the map Hi(V,Z) — ► H\(M C ,Z) induced by the inclu- 
sion V <^-> M c is surjective and consequently we can find 4 cycles v\, v 2 , ^3, ^4 
on the curve V such that 

A = ( L L "1 L 2 L 4 L 4 w i 

V L ^2 L W 2 / ^2 L U 2 



and La = {52t=i n k I j""* ^ j '■ n k £ Z}. The cycles ^1,^2,^3,^4 in P 
which we look for are 01,61,02,62 and they generate Hi(M c ,Z) such that 

A = 



/ai W l /fti^l /a2 W l 4 Wl 
/ ai W2 / 6l W2 J a2 W 2 J fe2 W2 
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is a Riemann matrix. We show that A = T* ,i.e., the period matrix of 
Prym*(V/£) dual of Prym(V/£). Consequently M c and Prym* (V / '£) 
are two abelian varieties analytically isomorphic to the same complex torus 
C 2 /La- By Chow's theorem [9], A c and Prym* (V / £) are then algebraically 
isomorphic. 

3.3 The Kowalewski rigid body motion 

The motion for the Kowalewski's top is governed by the equations 

to = to A Am, + 7 A Z, 7 = 7 A Ato, (39) 

where to, 7 and I denote respectively the angular momentum, the direc- 
tional cosine of the z-axis (fixed in space), the center of gravity which 
after some rescaling and normalization may be taken as / = (1,0,0) and 
Am = (toi/2, TO2/2, TO3/2) . The system (39) can be written 



TOi = TO 2 TO 3 , 

TO 2 = - TOl TO 3 + 273, 

TO 3 = -272, 



71 = 2 TO372 - TO273, 

7 2 = TO173 - 2to 3 7i, 

73 = m27l - TO472, 



(40) 



with constants of motion 
1 



Hi 



{m\ + m|) + to| + 271 = a, 

TO171 + TO272 + TO373 = c 2 , 



2 1 2 1 2 1 
= 7l + 72 + 73 = c 3 = 1, 



TOl + im 2 



(71 + ij2 




TOl — im 2 



(71 - 272) 



(41) 



c 4 . 



The system (40) admits two distinct families of Laurent series solutions : 



TOl (t) 



f + i(aj-2)a 2 + o (t) , 
- i (a\ - 2) a 2 + o (t) , 



t 



71 (t) 



(*), 



1 

w 



TO 2 (t) 



if - «f a 2 + o(t), 



72 (0 



2 

— ? 
2? 



(0, 



TO 3 (t) 



f + aia 2 + o (t) , 
=f + aia 2 + o (t) , 



73 (*) 



t+o(t), 



which depend on 5 free parameters ai, 05. By substituting these series in 
the constants of the motion Hi (41), one eliminates three parameters linearly, 
leading to algebraic relation between the two remaining parameters, which 
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is nothing but the equation of the divisor V along which the wij, 7$ blow up. 
Since the system (40) admits two families of Laurent solutions, then V is a 
set of two isomorphic curves of genus 3, V = T>\ + T>-\ : 

V £ : P (ai, a 2 ) = (a? - l) ((a? - l) - P (a 2 )) + c 4 = 0, (42) 

where P (a 2 ) = ci«2 — 2ec 2 a 2 — 1 and e = ±1. Each of the curve P e is a 
2 — 1 ramified cover («i,a 2 ,/3) of elliptic curves V £ : 

V° e : (3 2 = P 2 (a 2 ) - 4c 4 af, (43) 

ramified at the 4 points «i = covering the 4 roots of P(a 2 ) = 0. It was 
shown [12] that each divisor V £ is ample and defines a polarization (1,2) , 
whereas the divisor T>, of geometric genus 9, is very ample and defines a 
polarization (2,4). The affine surface M c = f] i=1 {Hi = C C 6 , defined 
by putting the four invariants (41) of the Kowalewski flow (40) equal to 
generic constants, is the affine part of an abelian surface M c with 

M c \ M c = V = one genus 9 curve consisting of two genus 3 
curves T> £ (42) intersecting in 4 points. Each 
V £ is a double cover of an elliptic curve T> £ (43) 
ramified at 4 points. 

Moreover, the Hamiltonian flows generated by the vector fields Xh 1 and Xh 4 
are straight lines on M c . The 8 functions 1 , f± = m\ , / 2 = m 2 , = 7713 , 

fi = 73 , h = fi + /I , h = 4/i/4 - 73/5 , h = (/271 - hi2)h + 2/472, 
form a basis of the vector space of meromorphic functions on M c with at 
worst a simple pole along T> Moreover, the map 

M c ~ C 2 / Lattice — > CP 7 , (h, t 2 ) » [(1, fi (h, t 2 ) , .... h (h, t 2 ))] , 

is an embedding of M c into CP 7 . Following the method (theorem 13), we 
obtain the following theorem : 

Theorem 14 The tori M c can be identified as M c = Prym*{V £ /V^, i.e., 
dual of Prym(T> £ /'D®) and the problem linearizes on this Prym variety. 



3.4 Kirchhoff's equations of motion of a solid in an ideal fluid 

The Kirchhoff's equations of motion of a solid in an ideal fluid have the form 
dH dH . dH dH , dH , dH 

Pi = P2^j — P3^r , h=P2x P3^ — m 2 ^ — h^-, 

dl 3 0I2 dp 3 dp 2 dl 3 0I2 

dH dH ■ dH dH dH dH , s 

P2 = P3~^ Pl^T , h=P3-R P17 h/ 3 ^ ilflT, (44) 

dh dls dpi dp 3 dh dls 

dH dH . dH dH dH , dH 

P3 = Pl-K-j P2-^T , h=Pl^ «2^~, 

di2 dli dp2 dpi dh oli 
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where (pi,P2,P3) is the velocity of a point fixed relatively to the solid, 
h, h) the angular velocity of the body expressed with regard to a frame 
of reference also fixed relatively to the solid and H is the hamiltonian. 
These equations can be regarded as the equations of the geodesies of the 
right-invariant metric on the group E (3) = SO (3) x R 3 of motions of 3- 
dimensional euclidean space M 3 , generated by rotations and translations. 
Hence the motion has the trivial coadjoint orbit invariants (p,p) and (p, l). 
As it turns out, this is a special case of a more general system of equations 
written as 

dH dH dH dH 

x = xA ^~ + v A ^~ ' y = y A + x A "s - ' 

ox oy ox oy 

where x = (xi, X2, x 3 ) £ M 3 et y = (2/1 ,2/2,2/3) £ R 3 - The first set can be 
obtained from the second by putting (x,y) = (l,p/e) and letting e — > 0. 
The latter set of equations is the geodesic flow on 50(4) for a left invariant 
metric defined by the quadratic form H. In Clebsch's case, equations (44) 
have the four invariants : 



Hi = H = - (a lP j + a 2 p\ + a 3 p\ + b x l\ + b 2 l\ + 6 3 /|) , 

H 2 = pl+pl+pl 
H 3 = 



Pih +P2I2+P3I3, 
Ha = \{bip\ + b 2 pl + b 3 pl + e {l\ + ll + ll)), 

with + S 2 ^- + = 0, and the constant g satisfies the conditions 

= 61(63-63) = 63(63-61) = b 3 ( bl -b 2 ) ^ The (44) can be written ■ th 

e 0,2 — tta (23— ai 01 — 0,2 J v ' 

form (11) with m = 6; to be precise 

dH 

x = f(x) = J— , x = {pi,p 2 ,p 3 ,h,l 2 ,k) J , (45) 






~P3 


P2 


P3 





-Pi 


~P2 


Pi 







where 



Consider points at infinity which are limit points of trajectories of the flow. 
In fact, there is a Laurent decomposition of such asymptotic solutions, 

x (t) = t' 1 (x(°> + x^t + x^t 2 + ...) , (46) 

which depend on dim(phase space) — 1 = 5 free parameters. Putting (46) 
into (45), solving inductively for the x^ k \ one finds at the th step a non- 
linear equation, x^ + f(x^) = 0, and at the k th step, a linear system of 
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equations, 

(L - kI)xW = { for fc = 1 

\ quadratic polynomial in x^ l \ for k > 1, 

where L denotes the jacobian map of the non-linear equation above. One 
parameter appear at the th step, i.e., in the resolution of the non-linear 
equation and the 4 remaining ones at the k th step, k = 1,...,4. Taking 
into account only solutions trajectories lying on the invariant surface M c = 

4 

P| {Hi (x) = Cj} C C 6 , we obtain one-parameter families which are param- 
eterized by a curve : 

V : 6 2 + ci/3 2 7 2 + c 2 aV + c 3 a 2 /3 2 + c 4 a/3 7 = 0, (47) 

where 6 is an arbitrary parameter and where a = x^ , (3 = x^ , 7 = 
Xg°^ parameterizes the elliptic curve 

£ : (3 2 = d\a 2 - 1, 7 2 = d\a 2 + 1, (48) 

with di, di such that: d\ + d\ + 1 = 0. The curve V is a 2-sheeted ramified 
covering of the elliptic curve £ . The branch points are defined by the 16 
zeroes of ci/3 2 7 2 + C2« 2 7 2 + csa 2 [3 2 + 040/^7 on £. The curve V is unramified 
at infinity and by Hurwitz's formula, the genus of V is 9. Upon putting 
( = a 2 , the curve V can also be seen as a 4— sheeted unramified covering of 
the following curve of genus 3 : 

C : (9 2 + Cl /3 2 7 2 + (c 27 2 + c 3 p 2 ) C) 2 - c 2 C/3 2 7 2 = 0. 

Moreover, the map r : C — ► C, (8, () (— 6>,£), is an involution on C and 
the quotient Co = C/r is an elliptic curve defined by 

C : V 2 = cj( (d 2 d 2 C 2 + (d 2 - d 2 ) C - 1) . 

The curve C is a double ramified covering of Co, C — ► Co, (6, 77, Q 1— > (r/, £), 

02 = _ Cl/3 2 7 2 _ ( C27 2 + C3/3 2) C + v 

v 2 = 4c {d\d 2 e + (di -d 2 )c-i). 

Let (ai, a-i, CL3, bi, 62, 63) be a canonical homology basis of C such that r (a±) = 
a 3, T (bi) = &3, t (0,2) = —a,2 and r(^) = — &2 for the involution r. Using 
the Poincare residu map, we show that 

d( (d( d( 

— — , — -— , u; 2 — 7-, 
77 6^77 6/77 

form a basis of holomorphic differentials on C and r* (wo) = ^o,r* (a;^) = 
— uJk (k = 1,2). The flow evolues on an abelian surface M c C CP 7 of period 

matrix ( ^ ^ ^ J , Im ( ^ ^ j > 0. Following the method (theorem 



C 



13), we obtain 
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Theorem 15 The abelian surface M c can be identified as Prym{C / C$) . 
More precisely 



4 



P| {x € C 6 , Hi (x) = q} = Prym{C/C )\V, 



i=i 



where V is a genus 9 curve (41), which is a ramified cover of an elliptic 
curve £ (48) with 16 branch points. 

4 Generalized algebraic completely integrable sys- 
tems 

Some others integrable systems appear as coverings of algebraic completely 
integrable systems. The manifolds invariant by the complex flows are cov- 
erings of abelian varieties and these systems are called algebraic completely 
integrable in the generalized sense. 

Consider the case F% = 0, (see section 3.1) and the following change of 
variables 

Zi = Qi, z 2 = q2, z 3 = p 2 , z 4 = piqi, z 5 = p\ - q\q\. 

Substituting this into the constants of motion Fi,Fi,Fs leads obviously to 
the relations 



whereas the last constant leads to an identity. Using the differential equa- 
tions (28) combined with the transformation above leads to the system of 
differential equations 



The last equation (28) for z$ leads to an identity. Thus, we obtain the po- 
tential constructed by Ramani, Dorozzi and Grammaticos [25,7]. Evidently, 
the functions H\ and H2 commute : {.Hi, #2} = 0. The system (50) is 
weight-homogeneous with qi,q2 having weight 1 and Pi,P2 weight 2, so that 
Hi and Hi have weight 4 and 8 respectively. When one examines all pos- 
sible singularities, one finds that it possible for the variable qi to contain 
square root terms of the type i 1//2 , which are strictly not allowed by the 
Painleve test (i.e., the general solutions have no movable singularities other 
than poles). However, these terms are trivially removed by introducing the 




(49) 




(50) 



qi (3qj + 8q 2 2 ) . 
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variables z±, . . . , z§ which restores the Painleve property to the system. Let 
B be the affine variety defined by 



B= f]{zeC 4 :H k (z)=b k }, (51) 
fc=i 

where (61,62) G C 2 . 

Theorem 16 a) The system (50) admits Laurent solutions in t 1 / 2 , depend- 
ing on 3 free parameters: u, v and w. These solutions restricted to the surface 
B(51) are parameterized by two copies Y\ and T_i of the same Riemann sur- 
face of genus 16. 

b) The invariant surface B(51) can be completed as a cyclic double cover B of 
the abelian surface A, ramified along the divisor C\+C-\. The system (50) is 
algebraic complete integrable in the generalized sense. Moreover, B is smooth 
except at the point lying over the singularity (of type A3) ofC\ +C_i and the 
resolution B of B is a surface of general type with invariants : X(B) = 1 
and p g (B) = 2. 

Proof, a) The system (50) possesses 3-dimensional family of Laurent solu- 
tions (principal balances) depending on three free parameters u, v and w. 
There are precisely two such families, labeled by e = ±1, and they are ex- 
plicitly given as follows 

11 5 13 7 

gi = -j=(u--uh + vt 2 -—u 7 t 3 + -u(ju 3 v-—u s + 5ew)t' 1 + •••), 

<72 = ^(^e-^eu 2 t + ^euH 2 + ^eu(^u 5 -3v)t 3 + wt 4 + •••), (52) 
pi = ^(--^ + 3^ 2 -^V + 

lu(^u 3 v-^u 8 + 3ew)t 4 + ■■■), 

P2 = ^(~e + ^eu A t 2 + ^eu(^u 5 -3v)t 3 + 3wt A + •••)• 

These formal series solutions are convergent as a consequence of the majorant 
method. By substituting these series in the constants of the motion Hi = b\ 
and H2 = 62, one eliminates the parameter w linearly, leading to an equation 
connecting the two remaining parameters u and v : 

f}JT AO Q 

r : ~r uv3 + TT^v 2 + TTTno (-9829u 8 + 26112fli) u 3 v (53) 
4 o4 8192 

10299 16 123 TT o TT 15362 98731 

According to Hurwitz' formula, this defines a Riemann surface V of genus 16. 
The Laurent solutions restricted to the surface -B(51) are thus parameterized 
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by two copies T_i and Ti of the same Riemann surface T. 
b) The morphism ip : B — ► A, (qi,q2,Pi,P2) 1 — ► {zi, z 2 , z 3 , Z4, z 5 ), maps 
the vector field (50) into an algebraic completely integrable system (1) in 
five unknowns and the affine variety -6(51) onto the affine part A (30) of an 
abelian variety A with A\A = C\ + C-\. Observe that p is an unramified 
cover. The Riemann surface T(53) play an important role in the construction 
of a compactification B of B. Let us denote by G a cyclic group of two 
elements {-1,1} on V£ = U£ x {r G C : < |r| < <5}, where r = t l/2 
and U£ is an affine chart of T £ for which the Laurent solutions (52) are 
defined. The action of G is defined by (—1) o (u,v,t) = (—u,—v,—t) and 
is without fixed points in V£ . So we can identify the quotient V£ /G with 
the image of the smooth map hi :V£^>B defined by the expansions (10). 
We have (—1, l).(u, v, r) = (— u, —v, r) and (1, —l).(u, v, r) = (u, v, — r), i.e., 
G xG acts separately on each coordinate. Thus, identifying V£ /G 2 with the 
image of p o hi in A. Note that Bl = V£ /G is smooth (except for a finite 
number of points) and the coherence of the Bl follows from the coherence of 
V£ and the action of G. Now by taking B and by gluing on various varieties 
-B|\{some points}, we obtain a smooth complex manifold B which is a double 
cover of the abelian variety A (constructed in proposition 2.3) ramified along 
C\ +C-i, and therefore can be completed to an algebraic cyclic cover of A. 
To see what happens to the missing points, we must investigate the image 
of T x {0} in UB£. The quotient T x {0}/G is birationally equivalent to the 
Riemann surface T of genus 7 : 

T . 65 3 93 3 2 J_ (_ 9829x 4 26112&i) x 2 y 
4 y 64 y 8192 v ' y 

( 10299 o 123, 4 , 15362 98731 \ 

where y = uv,x = u 2 . The Riemann surface T is birationally equivalent 
to C. The only points of T fixed under (u,v) 1— > (—u,—v) are the points 

2— 1 

at 00, which correspond to the ramification points of the map F x {0} — > 
T : (n, v) 1 ^ (x, y) and coincides with the points at 00 of the Riemann 
surface C. Then the variety B constructed above is birationally equivalent to 
the compactification B of the generic invariant surface B. So B is a cyclic 
double cover of the abelian surface A ramified along the divisor C\ + C_i, 
where C\ and C-\ have two points in commune at which they are tangent to 
each other. It follows that The system (8) is algebraic complete integrable 
in the generalized sense. Moreover, B is smooth except at the point lying 
over the singularity (of type .A3) of C\ + C-\. In term of an appropriate 
local holomorphic coordinate system (X,Y,Z), the local analytic equation 
about this singularity is X 4 + Y 2 + Z 2 = 0. Now, let B be the resolution 
of singularities of B, X{B) be the Euler characteristic of B and p g (B) the 
geometric genus of B. Then B is a surface of general type with invariants : 
X{B) = 1 and p g (B) = 2. This concludes the proof of the theorem. 
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Remark 4.1 The asymptotic solution (52) can be read off from (31) and 
the change of variable : q\ = yfz\,q2 = %2,Pi = £4/(71, f>2 = %3- The function 
z\ has a simple pole along the divisor C\ + C_i and a double zero along 
a Riemann surface of genus 7 defining a double cover of A ramified along 

Ci+C-l 
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